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ABSTRACT 

We compare the highly clustered populations of very high redshift galaxies with proto-clusters identified 
numerically in a standard ACDM universe (Slo = 0.3, Aq = 0.7) simulation. We evolve 256 3 dark matter par- 
ticles in a comoving box of side 150/i _1 Mpc. By the present day there are 63 cluster sized objects of mass in 
excess of 1O 14 /i _1 M in this box. We trace these clusters back to higher redshift finding that their progenitors at 
Z = 4-5 are extended regions of typically 20^4-0 Mpc (comoving) in size, with dark halos of mass in excess of 
IO'Vm© and are overdense by typically 1.3-13 times the cosmological mean density. Comparison with the 
observation of Lya emitting (LAEs) galaxies at z = 4.86 and at z = 4. 1 indicates that the observed excess clus- 
tering is consistent with that expected for a proto-cluster region if LAEs typically correspond to massive dark 
halos of more than 10 12 /i _1 M Q . We give a brief discussion on the relation between high redshift concentration 
of massive dark halos and present day rich clusters of galaxies. 
Subject headings: galaxies:clusters:general - cosmology:theory - methods:numerical 



1. INTRODUCTION 

The formation and evolution of structures in the universe, 
such as galaxy clusters and large-scale structures, are part 
of the most important issue in astrophysics. Since clus- 
tering properties of galaxies in the distant universe give 
us great clues to this issue, many observations have been 
done in orde r to detect (proto) clusters of galaxies at 
high redshi fts (Steid el et al.lll998Ll2000HCampos et al]|1999[ 
Pente ricci et alJl2000t IRhoads et alJl2000t lOuchi et alJboOll 
l2003riMiIevetaITl2004l) . These deep surveys of galaxies, 
which have significantly advanced our understanding of the 
properties and distribution of galaxies at high redshifts, are 
equally important in constraining the underlying structure for- 
mation scenarios. 

IVenemans et alJ ((2002) showed a region around a radio 
galaxy in which number density of Lya emitters (LAEs) is 
much higher than t he mean of the universe suggested by 
IRhoads et alJ (|2000) at z = 4.1. The region has size of 2.7 
x 1.8 Mpc (phys ical) an d mass of ~ 10 15 M Q . More recently, 
Shimasaku et alJ J2003> found that LAEs are clustered in an 
elongated region with a size of 20 x 50 Mpc (comoving) at 
z = 4.86, which is comparable to the size of present-day large- 
scale structures. In this elongated region, there is a circular 
region of high surface density of LAEs with 12 Mpc radius, 
which may be a progenitor of a galaxy cluster. They also esti- 
mated the bias parameter of the LAEs in the range of b ~ 3- 
16 for spatially flat low-density cosmological model (ACDM 
model). The elongated distribution of the LAEs in this region 
is also proposed to be a part of large-sca le structures, although 
observation by Shimasak u et al.l l|2£)04) of the same sky area 
which is closer to us by 3 9Mpc (comoving) does not show 
large structure. In addition, Havashino et al. (2004) reported 
that LAEs distribution at z = 3.1 in their observed area shows 
broad large scale structure. They claim that this structure can- 
not be explained in the context of the standard CDM model. 
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lOuchi et"aT] ( 1200 5) observed LAEs at z = 5.7 and discovered 
filamentary structures and voids. 

Although these regions which have much high number den- 
sities of LAEs are claimed to be proto-clusters, there is not 
yet enough detailed study on these possibilities. Therefore, 
we investigate formation and evolution of galaxy clusters in 
early universe using cosmological simulation. In this let- 
ter we report some properties of simulated proto-clusters at 
z = 5 and compare our n um erical results with obse rvations by 
IShimasaku et alJ J2003I) and lVenemans et alJ J2002I) . 

2. METHOD 

2.1. Numerical Method 

We perform a cosmological N-body simula- 
tion with Particle-Partic le-Particle-Mesh (P 3 M) 
algorithmf Hocknev & Eastwoodll9 81[). Our simulatio n code 
is the same one which was used in lYoshikawa. Jing. & Sutol 
(2000). We use the following values as parameters of our 
simulation: the Hubble constant in units of 100km s" 1 Mpc" 1 , 
h = 0.7, the density parameter, flo = 0.3, the baryon density 
parameter, £1/, = 0.015/i~ 2 , the root mean square density 
fluctuation amplitude on a scale 8/i -I Mpc, erg = 1.0, the 
power-law index of the primordial density fluctuation, 
n = 1.0, and the cosmological constant parameter, Ao = 0.7. 
We employ A^dm = 256 3 dark matter particles and the mass 
of each one is 2.15 x 1O 1O M . The size of the comoving 
simulation box, L\, ox , is 150/i _1 Mpc, and the box is on the 
periodic boundary condition. This size of the box is large 
enough to realize a sufficient number of clusters and large- 
scale structures at z = 0. We use the spline (S2) softening 
function for gravitational softening, and the softening length, 

e grav , is set to be L^ ox /(l0N^) (~ 60/i _1 kpc in the comoving 
scale). 

2.2. Identification of Dark Halos and Proto-clusters 

We identify da rk halos in a manner similar to that of 
ISuwa et all (120031) . A brief explanation of the identification 
is as follows: First, we define densities of dark matter par- 
ticles using an interpolation technique in the sa me way of 
the S moothed Particle Hydrodynamics method (Monaghan 
119921) . and pick up particles whose densities are more than 
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the virial density. Next, we perform the hierarchical friends- 
of-friends (HFoF) method ( Klvpin et al. 1999) for the selected 
dense particles. The maximum linking length Z max in HFoF 
method is defined as 0.21, where / is the mean inter-particle 
distance for all (not only dense) particles. Then we draw a 
sphere of which the center is on the densest particle of the 
group and seek the radius in which the mean density of to- 
tal matter is equal to the virial density. We regard the set 
of particles in the sphere as a virialized object. If a centroid 
of a sphere exists in another sphere, only the set which be- 
longs to the more massive sphere is used. We regard dark 
halos which have mass of more than 10 i4 /i _1 Mq at z = as 
galaxy clusters and we analyze properties of their progeni- 
tors to compare with the observations of proto-cluster regions 
dVenemans et alJ20 02: Shimasak u et alJ2003l) . 

We seek proto-clusters in high redshift universe as follows: 
First, we pick up all particles in each cluster at z = 0. Next, 
we trace those particles back to high-z epoch, e.g. z = 5. Fi- 
nally, we set a minimum cubic region which covers all of the 
particles. We call the region as "proto-cluster region" for the 
cluster. This region contains the proto-cluster of the present 
cluster. 

In order to investigate galaxy distribution in the proto- 
cluster regions, we identify galaxy-size dark halos (M > 
1O 12 /i~'M0) in the regions at z = 5. Dark halos of this mass 
scale are suggested to conta in LAEs and LBGs f rom their spa- 
tial distribution arguments ( Haman a"et all2004T) . 

2.3. Overdensity 

In proto-cluster regions, it is expected that dark halos which 
correspond to galaxies concentrate much denser than the 
background. Therefore the excess of number density of dark 
halos, <5haio, should be useful quantities to compare the result 
of the simulation to observations: 



"halo.BG 

where «haio is the halo number density averaged in the suit- 
able scale (e.g. 25 Mpc) and «haio,BG is the background halo 
number density. 
The mass overdensity, 6 mass , is defined in a similar way: 

^mass = 1 : (2) 

Pbg 

where p is the density of the region under consideration 
and peg is that of the background. The combination of 
the halo overdensity and the mass overdensity give us im- 
portant information of mass concentration, so-called bias 
information, which has been i nvestigated in analytic al and 
numerical works (e.g . iKaiserl 11984 iTaruva & Sutol l2000t 
lYoshikawaetall2001l) . 

3. RESULTS 

We find 63 galaxy clusters with M > 10 14 !i~ 1 Mq in the sim- 
ulation box at the present epoch (z = 0). We seek their progen- 
itors, i.e. proto-clusters, in high redshift universe and study 
their properties. We show an example of present galaxy clus- 
ters andproto-cluster regions in our simulation box at z = 5 in 
Figure^ 

In our numerical results, the size of the proto-cluster re- 
gions is in the range of 20^-OMpc in the comoving scale. 
The exten sion is similar to that of the LAEs dense region ob- 
served by Shimasaku et al. (2003). Several dark halos with 
M > 10 12 /! _I M Q already exist in proto-cluster regions at z = 5 




FIG. 1 . — Distributions of dark matter and dark halos in an example cluster 
and the corresponding proto-cluster region in the simulation box. Panel (a) 
and (b) show the distribution of dark matter at z = and z = 5, respectively. 
Panel (c) shows the distribution of dark halos which have mass of more than 
10 12 /! -1 Mq in the proto-cluster region shown in the panel (b). x- and y-axes 
are comoving position coordinate in the simulation. 
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FIG. 2. — Histogram of halo number overdensities, <5h a i , °f the proto- 
cluster regions (red solid line) and the randomly selected region (green 
dashed line) at z= 5, 



as shown in Figure^c). We obtain bulk velocity of each dark 
halo in the proto-cluster regions and found that dispersion of 
those are typically ~ 200km s" 1 . 

We calculate <5haio and <5 mass for each proto-cluster regions at 
z = 5, by assuming a smoothing scale of 25 Mpc (comoving) 
that is a typical proto-cluster size in our numerical results. In 
order to show that the values of these indicators have signifi- 
cant excess than the mean value of the universe, we randomly 
select (25Mpc) 3 regions in the simulation box and calculate 
<5haio and (5 mass for these regions. The number of the randomly 
selected regions is 630, which is 10 times more numerous than 
the proto-cluster regions. We show the histograms of cWo in 
Figure |3 where the red solid line is for the proto-cluster re- 
gions and the green dashed line is for the randomly selected 
regions. The red solid line is significantly different from the 
green dashed line. Figure [2] shows that the typical value of 
(Wo is ~ 3 for the proto-cluster regions, and its variance is 
very large. On the other hand, Jhaio in most of the randomly 
selected regions are almost -1, where this value means that 
there is no massive dark halo in the regions. 

The mass overdensities, S mass , for the proto-cluster regions 
are clearly different from the randomly selected regions as 
shown in Figure [3] <5 mass of the proto-cluster regions are in 
the range of 0.2-0.4, while those of the randomly selected 
regions range from -0.2 to 0.2. The typical value of i5 mass of 
the proto-cluster regions is ~ 0.4 and the variance of 5 mass is 
much smaller than that of (5h a io- 

In order to show the bias of dark matter halo distribution in 
the proto-cluster regions, we plot <5 mass and <5h a io for the proto- 
cluster regions and the randomly selected regions in Figure 
|4] Red triangles, blue crosses, and green squares indicate the 
proto-cluster regions, the randomly selected regions which 
overlap more than 50% with some proto-cluster region, and 
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FIG. 3. — Same as Figure|2] but for 5 m ass- 




FIG. 4. — Correlation map of (5 raaS s and <5[, ; ,i . Red triangles, blue crosses, 
and green squares indicate the proto-cluster regions, the randomly selected 
regions which overlap more than 50% with some proto-cluster region, and 
other randomly selected regions, respectively. The analytical relation is in- 
dicated by a thick black line. Green and red lines are correspond to average 
value of simulated regions randomly selected (both field and overlapped) and 
proto-clusters, respectively. 



other randomly selected regions, respectively. Figure|4]shows 
correlation between the values of <$haio and S mass for c5 mass > 0, 
although the dispersion is very large. It is clear that the re- 
gions with large <5 mass have large <5haio- This relation corre- 
sponds to the bias of dark matter halo distribution in the proto- 
cluster regions. 

We compare the <5 mass -t5haio relation in our numerical re- 
sults wifhjh^jmalyjical result given by the natural bias the- 
ory iMo & Whitd[l9 96) in Figure |4] In Figure we show 
the analytical relation by a thick black line, where we use 
the Press-Schchter mass function formula and the relation of 
linear extrapol ation to nonlinear evolu tion for density pertur- 
bation dCarroll. Press. & Turnerin"992l) . Green and red lines 
in Figure 0] indicate average value of the randomly selected 
regions and the proto-clusters, respectively. The relation be- 
tween the t> mass -<5haio given by the natural bias theory agrees 
well with our numerical results as shown in Figure |4] ex- 
cept for gradual deviation <5 mass > 0.4. This deviation may 
be explained by the effect of limitation of our simulation box 
size. Our simulation box size is not enough to consider large- 
scale (lOO/r'Mpc) component of density fluctuation and this 
limitation may cause underestima tion of mass function for 
10 12 M Q dark halos at z = 5 JBagla & Ravl2005l) . 

It is very important to show a critical value of 5haio> <5haio,c, 



CO 

d 



CD 

d 




2 4 6 8 10 



FIG. 5. — Fraction of randomly selected cubic regions which contain a 
galaxy cluster at z = as a functions of <5|, a i at z = 5. The size of the cubic 
regions is 25 Mpc (comoving). 



by which we can select a proto-cluster region. In order to 
search <5haio,c, we compare values of halo overdensity, (5h a io, 
in the randomly selected regions at z = 5 and masses of the 
largest dark halos in the regions at the present epoch. Figure|5] 
shows a fraction of the regions which contain galaxy clusters 
with M > lO'VM© at z = as a function of (5haio at z = 5. 
It is clear in Figure |5]that most regions with Jhaio > 3 at z = 5 
will evolve rich clusters with mass of more than 10 l4 /i _1 M© 
at z = 0. 

4. DISCUSSION 

We have investigated proto-cluster regions and the differ- 
ence between those regions and background field regions us- 
ing N-body cosmological code. In order to investigate proto- 
clusters, we identify 63 galaxy clusters at in our cosmological 
simulation z = 0, and trace dark matter particles which belong 
to those, back to a high redshift epoch. We study properties 
of the proto-cluster regions by analyzing the mass overden- 
sity, <5 maS s, and the halo overdensity, <5h a i , °f which definitions 
are given in ^2] The results at z = 5 are shown in FiguresEIEl 
and|4] As expected from the bias theory (e.g. lTaruva & Suto 
2000), <5h a io shows excess over <5 mass and variance of Jhaio is 
much larger than that of <5 mass . 

In order to show difference between the proto-cluster re- 
gions and the field, we calculate (5 mass and <5h a i for the ran- 
domly selected regions with the same size of the typical proto- 
cluster region at z = 5 (25 Mpc in the comoving scale). We 
select 630 regions (10 times greater than the proto-cluster re- 
gions) for the randomly selected regions. Clustering proper- 
ties of the dark halos in these regions are significantly differ- 
ent from those in the proto-cluster regions. The halo over- 
densities for most of the randomly selected regions are much 
smaller than those for the proto-cluster regions as shown in 
Figure |3 

It is very interesting to study a critical value of i5haio by 
which we can select a region as the proto-cluster regions. 
From Figure [5] we have found that more than 80% of the 
regions with i5h a io > 3 at z = 5 contain rich clusters (M > 
10 14 /! _1 Mq) at the present epoch. Thus, we conclude that 
^haio > 3 is a good criterion to distinguish proto-cluster re- 
gions from field regions at z = 5. This criterion is very useful 
to judge whether observed galaxies excess regions are actual 
proto-clusters or not. 

It should be pointed out that there are a few regions which 
do not contain clusters at z = but have large i5haio values 
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(3-7) at z = 5, although a fraction of those regions are very 
small. We discuss that the physical origin of the large 5haio 
is non-linearity and stoch asticity of bias parameters proposed 
by Taruva & Suto (2000). They studied variance of biasing 
parameters, and conclude that the variance increases strongly 
with redshift. They also conclude that stochasticity of the bi- 
asing is generate d by the scatter in the halo mass distribution 
at higher redshift. Yoshik awa et alJ ( 120011) confirmed these re- 
sults for < z < 3 by large P 3 M simulation. We conclude that 
our results for z = 3 are consistent with their results and the 
trend of evolution of bias variance is still true for 3 < z < 5. 

We compare our results with the observation of 
Shimasaku et alJ (|2003). Properties of the LAEs con- 
centrated region reported by them are as follows: (1) 
Diameter of the region is 25 Mpc (comoving unit), (2) 
Projected overdensity of LAEs 5^ ~ 2, (3) Bias parameter 
is estimated 3-16 for ACDM model, and the best-fit value 
is b ~ 6, i.e. <5 mass ~ 0.3. (4) The number of LAEs in the 
region is about 20. Properties from (1) to (3) are consistent 
with our numerical results of proto-cluster regions, if LAEs 
correspond to dark h alos with mass more than IO^/t'Mq, 
which is suggested by Hamana et al. (2004) using correlation 
function on small-scales. The typical number of dark halos in 
the simulated proto-clusters is 5-10, which is about half of 
that of observed LAEs. One explanation of this discrepancy 
is due to possibility that some dark halos have more than one 



LAEs. We speculate t hat some pairs of LAEs in Fig. 3 of 
Shimasaku et al. (2003) are included in the same dark halo. It 
is also possible that some of observed sampl e of LAE s in the 
region are low redshift interlopers. Shima saku et alJ ( 120031) 
e stimate the contaminatio n of their sample to be about 20%. 

Venemans et al. (2002) also reported properties of a LAEs 
rich region. The size of the region is ~ 14 x 10 Mpc (comov- 
ing unit). They estimate the region is overdense i n LAEs by 
a fact or of 15 compared with the blank field (Rhoads et al. 
2000). The overdensity is much larger than that of proto- 
clusters at z = 4 in our simulation and the size of their observed 
region is smaller than the typical size of the proto-clusters in 
our simulation. Therefore, we suggest that the region ob- 
served by them is not a whole proto-cluster region, but the 
central region of the proto-cluster because of small size of the 
observed region and their high value of the overdensity. 
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